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RESEARCH  SUMMARY 

Distributions  of  Douglas-fir  beetles  and  entomophagous 
insects  were  sampled  weekly  at  3  m  intervals  on  stems  of 
trees.  Regression  models  were  developed  for  predicting 
abundance  of  Douglas-fir  beetles  in  relation  to  entomo- 
phagous species,  date,  and  sample  characteristics.  The 
logistic  function  was  used  to  define  the  probability  of  the 
presence  of  entomophagous  insects  on  samples.  An  aver- 
age of  58  percent  mortality  of  Douglas-fir  beetle  progeny 
was  attributed  to  entomophagous  insects.  Suggestions  are 
given  on  locations  of  samples  and  number  of  samples  need- 
ed for  a  given  precision  of  estimate  for  Douglas-fir  beetle 
and  entomophagous  species. 
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INTRODUCTION 

Infestations  of  the  Douglas-fir  beetle  Dendroctonus 
pseudotsugae  Hopk.  cause  great  economic  losses  of  its 
principal  host,  Pseudotsuga  menziesii  [Mirb.]  Franco 
(Furniss  and  Orr  1978).  Mature  trees  are  commonly  killed 
after  predisposing  factors  such  as  windstorms  and 
droughts  enable  the  beetle  population  to  increase. 

The  Douglas-fir  beetle  develops  through  its  four  life 
stages  (fig.  1)  hidden  under  the  bark  of  its  host  tree.  Most 
beetles  overwinter  as  callow  adults,  then  fly  and  infest  other 
Douglas-fir  trees  during  spring.  The  females  bore  into  the 
bark  and  excavate  unbranched  egg  galleries  upward  in  the 
phloem  for  an  average  distance  of  1 5  to  20  cm.  The  density 
of  egg  galleries  ranges  from  5  to  12  per  0.1  m^  in  standing 
trees.  Eggs  are  laid  alternatively  on  opposite  sides  of  the 
gallery  at  a  rate  of  3.1  per  cm  in  Idaho.  Eggs  hatch  in  a  week 
or  two.  The  larvae  then  mine  fanlike  (fig.  2)  from  the  egg 
gallery.  Larvae  complete  their  mines  in  a  month  or  so,  then 
construct  cells  in  which  they  transform  to  pupae  for  a  brief 
time  before  becoming  adults.  Death  of  thetree resultsfrom 
a  combination  of  girdling  by  larvae  and  the  action  of  micro- 
organisms associated  with  the  beetle. 

Distribution  of  beetles  and  entomophagous  insects  in 
trees  was  investigated  in  central  Idaho  (Kline  1 963)  as  a  part 
of  a  long-term  study  of  population  dynamics  by  Furniss. 
That  thesis  presented  graphically  the  abundance  of  these 
organisms  by  sampling  height  and  date.  In  this  paper,  we 
utilize  regression  analyses  and  the  logistic  function  to  re- 
late quantitatively  the  Douglas-fir  beetle  population  density 
to  the  number  of  predators  and  parasites  in  relation  to  time 
and  location  in  trees.  The  improved  methods  provide  a 
means  of  expressing  relationships  that  previously  were 
either  not  recognized  or  incapable  of  definition.  This  infor- 


mation will  facilitate  the  evaluation  of  Douglas-fir  beetle 
infestations  for  such  purposes  as  predicting  trend  of  popu- 
lation and  damage,  relating  Douglas-fir  beetle  populations 
to  factors  affecting  tree  susceptibility;  and  will  aid  in  de- 
veloping improved  sampling  procedures. 


Figure  1  .--Douglas-fir  beetle  life  stages. 
(A)  adult;  (B)  eggs  in  gallery;  (C) 
larva;  (D)  pupa. 
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Figure  2. --Douglas-fir  beetle  egg  gallery 
and  larval  mines. 


METHODS 

The  study  was  scheduled  to  begin  when  sufficient  time 
had  elapsed  after  beetle  attack  to  reveal  trees  that  were  go- 
ing to  die.  We  thereby  eliminated  tree  resistance  as  a 
source  of  mortality  to  beetle  progeny,  and  galleries  were 
nearly  completed,  which  resulted  in  less  variation  in  weekly 
measurements. 

Five  infested  trees  were  felled  weekly  for  sampling  begin- 
ning June  30, 1960.  Sampling  ended  September  2.  Samples 
(fig.  3)  were  taken  at  3  m  intervals,  beginning  1.5  m  above 
ground  to  a  diameter  of  20  cm.  Numbers  of  entomo- 
phagous  insects  were  measured  on  a  15cm  x  30cm  sample, 
but  because  Douglas-fir  beetle  galleries,  ventilation  holes, 
and  progeny  were  more  abundant,  they  were  measured  on 
an  adjacent  93  cm^  circular  sample  cut  with  a  bark  punch 
(Furniss  1962a).  Other  data  obtained  at  each  sample  loca- 
tion were:  diameter  of  trunk,  bark  thickness,  and  height  of 
sample. 


Figure  3. —Comparison  of  samples. 
Circular  sample  was  used  to  measure 
Douglas-fir  beetle  galleries  and 
progeny;  rectangular  sample  was 
used  to  count  entomophagous  insects 
on  adjacent  area. 

Analyses 

We  used  actual  sample  values  ratherthan  converting  data 
to  a  common  area  basis.  Thus,  all  counts  of  entomo- 
phagous insects  are  on  a  15  cm  x  30  cm  (0.45  m^)  basis 
whereas  densities  of  Douglas-fir  beetle  progeny  and  length 
of  egg  galleries  are  on  a  93  cm^  basis. 

To  estimate  rates  of  survival  of  Douglas-fir  beetles  from 
predation  and  parasitism,  we  used  ratios  of  numbers  of 
Douglas-fir  beetle  progeny  to  the  total  length  of  egg  gal- 
leries. Initial  beetle  population  (eggs)  is  related  to  egg  gal- 
lery length  (Furniss  1957),  there  being  an  average  of  3.1 
eggs/cm. 

Simple  linear  regression  was  used  to  relate  Douglas-fir 
beetle  population  indices  to  entomophagous  insects.  To 
determine  the  variation  explained  by  regression  models,  we 
calculated  multiple  regression  coefficients  (R^). 
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The  following  list  of  variables  was  considered  a  basic  set 
for  the  presence  or  absence  of  entomophagous  insects: 
date,  height,  length  of  egg  galleries,  maximum  bark  thick- 
ness, diameter,  and  number  of  ventilation  holes.  Trans- 
formations of  those  variables  were  also  tested  in  the  model 
building. 

In  developing  the  model,  the  response  was  fitted  to  the 
entire  list  of  variables  first.  Coefficients  in  the  model  were 
tested  for  significance  by  Student's  t-test.  Variables  with 
nonsignificant  (a  <  0.05)  coefficients  were  deleted  and  the 
coefficients  for  the  remaining  variables  were  recomputed. 

The  logisticfunction  (Walkerand  Duncan  1967)  was  used 
to  define  the  probability  of  the  presence  of  entomophagous 
insects  using  Hamilton's  (1 974)  program.  The  goodness  of 
fit  of  the  predicted  and  observed  values  was  tested  by  chi- 
square  analysis. 

Relationships  of  dependent  (e.g.,  number  of  Douglas-fir 
beetles)  and  independent  (e.g.,  number  of  entomophagous 
insects)  variables  were  portrayed  with  an  HP9820  calcula- 
tor/plotter. The  same  device  was  used  to  plot  the  logistic 
function  relating  entomophagous  insects  {Coeloides  or 
Medetera)  and  sample  characteristics. 

The  efficiency  of  sampling  beetle  progeny  at  1.5  m  and 
4.6  m  above  ground  was  compared  by  the  degree  to  which 
their  means  were  correlated  with  tree  means.  Numbers  of 
samples  needed  to  estimate  average  gallery  length,  beetle 
progeny,  or  entomophagous  insects  were  calculated  for  a 
0.20  coefficient  of  variation  of  the  mean  (CVj^)  using  the 
formula  n  =  (CVx/0.20)2.  This  is  equivalent  to  setting 
sample  size  for  a  desired  relative  variance  (Kish  1965). 


BIONOMICS  OF  SPECIFIC  PREDATORS 
AND  PARASITES 

The  insect  enemies  of  the  Douglas-fir  beetle  found  in  this 
study  are  discussed  in  the  general  order  of  their  density  on 
samples. 

COELOIDES  VANCOUVERENSIS  [Dalla  Torre] 

In  thin  bark  portions  of  trees,  C.  vancouverensis 
{=brunneri  Viereck)  can  be  very  effective  in  parasitizing  a 
large  percentage  of  the  Douglas-fir  brood  (fig.  4B)  (Ryan 
and  Rudinsky  1962).  For  instance,  Bedard  (1933)  reported 
an  average  of  29  percent  of  beetle  larvae  were  parasitized, 
varying  from  0-10  percent  in  the  base  to  70-95  percent  in  the 
top. 

Females  deposited  eggs  through  the  bark  (fig.  4A)  onto 
beetle  larvae,  beginning  in  early  July.  By  theend  of  July  the 
parasite  was  found  frequently  in  trees.  Due  to  a  different 
rate  of  development,  some  parasite  progeny  overwinter  as 
larvae.  Others  mature  and  emerge  in  August  to  lay  eggs  on 
other  beetle  larvae.  The  proportion  of  fast-developing 
adults  may  be  95  percent  in  western  Oregon  (Ryan  and 
Rudinsky  1962)  or  a  small  minority  in  northeastern  Wash- 
ington (Bedard  1933).  Each  C.  vancouverensis  larva  con- 
sumes only  a  single  beetle  larva,  so  the  numiber  of 
C.  vancouverensis  larvae  (or  their  cocoons  which  are  easily 
counted)  equals  the  number  of  bark  beetle  larvae  that  have 
been  parasitized. 


Figure  4.— Coeloides  vancouverensis 
ovipositing  (A)  on  beetle  larva  in  bark, 
and  cocoons  (B)  marking  demise  of 
beetle  larvae. 


MEDETERA  SPP. 

Both  Medetera  aldrichii  Wheeler  and  M.  vidua  Wh.  were 
present,  but  the  latter  was  much  less  abundant.  Adults 
(fig.  5A)  began  to  appear  in  late  May,  becoming  abundant  in 
June.  Eggs  are  laid  on  small  groups  of  1  to  3  in  bark  crevices 
(Bedard  1933).  Larvae  (fig.  5B)  began  to  appearon  samples 
June  30.  The  average  number  of  Douglas-fir  beetle  larvae 
killed  per  Medetera  larva  is  unknown. 
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Figure  5.-Medetera  adult  (A)  and 
larva  (B). 


ENOCLERUS  SPHEGEUS  Fab. 

Enoclerus  sphegeus  Fab.  adults  (fig.  6A)  emerged  in  May 
and  preyed  on  attacking  adult  Douglas-fir  beetles.  Eggs  are 
laid  in  clusters  under  scales  of  the  outer  bark  (Kline  and 
Rudinsky  1964).  Thenewly  hatched  larvae  apparently  enter 
the  beetle  galleries  through  holes  made  by  the  bark  beetles, 
after  which  they  feed  on  immature  stages  of  D. 
pseudotsugae. 

Larvae  of  E.  sphegeus  (fig.  6B)  were  present  on  samples 
throughout  the  sampling  period  (fig.  7)  although  they 
diminished  in  abundance  after  mid-July  when  most  larvae 
migrated  to  the  root  crown  where  they  formed  pupal  cells  in 
the  outer  bark  or  duff.  In  laboratory  tests,  larvae  ate  an 
average  of  0.5  D.  pseudotsugae  larva  per  day  (Bedard  1 933) 
and  15  to  38  Dendroctonus  ponderosae  Hopk  larvae,  de- 
pending on  prey  size  (Amman  1970).  Such  high  mortality 
rates  seem  unlikely  in  nature  due  to  the  somewhat  dis- 
persed and  isolated  nature  of  prey  larvae. 
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Figure  6.— Two  predators  of  similar  habit:  Enoclerus  sphegeus  adult  (A) 
and  larva  (B);  and  Temnochila  chlorodia  adult  (C)  and  larva  (D)  search- 
ing for  prey. 
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Figure  7. --Average  number  of  Enoclerus  sphegeus  per  sample  per  tree  by 
day  of  the  year  on  which  the  tree  sampled. 


THANASIMUS  UNDATULUS  Say 

Larvae  of  Thanasimus  undatulus  Say  were  present  at  low 
frequency  (fig.  8).  They  may  prefer  smaller  bark  beetles 
such  as  Pseudohylesinus  nebulosus  (Lec.)  or  Scolytus 
monticolae  (Sw.)  (Schmitz  1978).  The  seasonal  history  is 
similar  to  that  of  E.  sphegeus  with  the  exception  that  the 
larvae  do  not  migrate  to  the  base  of  the  tree. 

TEMNOCHILA  CHLORODIA  [Mann.] 

T  chlorodia  larvae  (fig.  6D)  were  found  throughout  the 
summer  in  abundance  similar  to  the  clerids  (fig. 9).  More 
recent  work  indicates  that  this  metallic  green  or  blue  preda- 
tor (fig.  6C)  is  much  less  abundant  in  northern  Idaho 
(Furniss  and  others  1979).  It  also  seems  less  important  in 
Douglas-fir  of  western  Oregon  (Kline  and  Rudinsky  1964). 
The  seasonal  history  of  T.  chlorodia  follows  that  of  E. 
sphegeus  except  that  the  larvae  do  not  migrate. 


ROPTROCERUS  ECCOPTOGASTERI  Ratz, 
CECIDOSTIBA  BURKEI  Crawford,  and  CECIDOSTIBA 
DENDROCTONI  Ashmead 

These  pteromalids  are  grouped  together  because  the 
larvae  (fig.  10)  are  similar  in  appearance  and  habit.  They 
may  be  more  important  as  parasites  of  Pseudohylesinus 
nebulosus  [Lec]  and  Scolytus  monticolae  [Sw.j  in  tops  of 
trees  killed  by  D.  pseudotsugae.  These  parasites  were  most 
abundant  in  the  late  summer  (fig.  11). 

BELOSTA  ALBIPILOSA  Hardy 

The  scenopinid  fly,  B.  albipilosa  [Hardy]  (fig.  12)  was  dis- 
covered in  this  study  for  the  first  time  preying  on  Douglas- 
fir  beetle  progeny  beneath  bark  (Kelsey  1969,  p.  275).  The 
segmented  white  larvae  are  elongate  and  possess  a  beak- 
like extension  of  the  head  (fig.  12).  They  were  recorded  on 
samples  only  in  late  summer  (fig.  13). 
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Figure  8. --Average  number  of  Thanasimus  undatulus  per  sample  per  tree 
by  day  of  the  year  on  which  the  tree  was  sampled. 
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Figure  9. —Average  number  of  Temnochila  chlorodia  per  sample  per  tree  by 
day  of  the  year  on  which  the  tree  was  sampled. 


Figure  ^0  .--Roptrocerus  larva 
recognizable  by  crescent  shape  and 
absence  of  cocoon. 
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Figure  11. —Average  number  of  Roptrocerus  and  Cecidostiba  per  sample  per 
tree  by  day  of  the  year  on  which  the  tree  was  sampled. 


Figure  ^2.—Belosta  albipilosa  adult  (A) 
and  larva  (B).  The  snakelike  larva  was 
discovered  as  a  predator  of  Douglas- 
fir  beetle  during  this  study. 
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Figure  1 3. --Average  number  of  Belosta  albipilosa  per  sample  per  tree  by  day  of 
the  year  on  which  the  tree  was  sampled. 
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RESULTS  AND  DISCUSSION 

After  June  20,  the  length  of  egg  galleries  pertree  did  not 
increase.  Therefore,  if  it  were  not  for  mortality  factors  the 
beetle  population  would  have  remained  constant  through 
the  summer;  however,  density  of  beetle  progeny  per  sample 
declined  (fig.  14). 

We  hypothesize  that  the  principal  cause  of  reduction  in 
beetle  numbers  in  a  successfully  infested  tree  was  due  to 
entomophagous  insects,  not  tree  resistance  or  competition 
among  progeny.  This  premise  issupported  by  thefollowing 
facts.  We  sampled  only  successfully  attacked  trees,  that  is, 
those  low  in  resistance.  Mortality  from  predation  and  para- 
sitism acted  to  thin  out  progeny,  lessening  competition. 
The   presence   of  the  braconid  wasp,  Coeloides 


vancouverensis  [Dalla  Torre]  (=  brunneri  Vier.)  coincided 
with  a  sudden  reduction  in  the  beetle  population.  The 
average  number  of  Douglas-fir  beetle  progeny  was  4.72  per 
decimenter  of  egg  gallery  per  tree.  For  sample  trees  with 
low  numbers  of  entomophagous  insects  the  average  number 
of  beetle  progeny  was  11 .26  per  dm  of  egg  gallery.  The  ratio 
was  32.48  for  the  only  tree  lacl<ing  entomophagous  insects 
when  sampled. 
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Figure  14.— Average  number  of  Douglas-fir  beetles  per  sample  per  tree  by  day 
of  the  year  on  which  the  tree  was  sampled. 


9 


Although  the  length  of  egg  galleries  did  not  increase  on 
samples  through  the  summer,  it  varied  greatly  from  tree  to 
tree.  To  remove  most  of  this  variation,  we  plotted  over  time 
the  number  of  beetles  per  average  length  of  egg  gallery 
(fig.  15). 

As  expected,  the  number  of  Douglas-fir  beetle  progeny 


were  inversely  related  to  numbers  of  entomophagous  in- 
sects (fig.  16).  Density  of  beetle  progeny  at  low  numbers  of 
any  observed  enemies  was  1 1.26/dm  of  egg  gallery  but  only 
2.76/dm  (76  percent  fewer)  at  high  frequency  of  parasites 
and  predators.  Progeny  averaged  4.72/dm  (58  percent 
fewer). 
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Figure  15. --Average  number  of  Douglas-fir  beetles  per  dm  of  egg  gallery  per 
^  tree  by  day  of  the  year  on  which  the  tree  was  sampled. 
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Figure  16. —Average  number  of  Douglas-fir  beetles  per  dm  of  egg 
gallery  per  tree  vs.  average  number  of  entomophagous  insects 
per  sample  per  tree. 
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Population  Models 

DOUGLAS-FIR  BEETLES 

As  a  beginning  point  we  used  simple  linear  regression 
analysis  to  answer  the  question:  Is  the  length  of  egg  gallery 
per  sample  a  function  of  tree  diameter,  sample  height,  and 
maximum  bark  thickness? 

Five  trees  were  selected  at  random  from  each  of  five  dia- 
meter breast  height  classes  (table  1).  A  quadratic  effect  of 
maximum  bark  thickness  is  suggested  by  these  data,  the 
average  length  of  egg  galleries  being  greatest  at  a  maxi- 
mum bark  thickness  of  1.27  to  1.78  cm.  A  quadratic  form  is 
also  used  for  the  sample  height  in  the  regression  model: 

where: 


Y^       =  length  of  egg  gallery  in  the  /th  sample, 

xi  =  diameter  at  1.5  m  oftreefrom  which /'th  sample  is 

'  taken, 

X2  =  maximum  bark  thickness  for  the  /th  sample, 

X3j  =  height  at  which /th  sample  occurred  on  the  tree, 

e,       -  error  associated  with  /'th  sample. 

This  regression  model  was  then  fitted  to  the  total  data  set. 

The  regression  model  accounted  for  only  a  small  per- 
centage of  the  variation  in  the  average  length  of  egg  gal- 
leries as  indicated  by  the  multiple  correlation  coefficient, 
=  0.111.  The  regression  coefficient/significance  levels 
were: 


bo  =  6.685/0.001  bA  =  0.081 8/.  054 

bl  =  -0.1109/.015  ^5  =  0.008/.  047 

t)2  =  -1.192/0.485 

^3  =  0.2691/0.627 

The  significant  coefficients  in  the  model,  i.e.,  those 
statistically  different  from  zero,  are  tree  diameter  and 
sample  height.  The  coefficients  for  maximum  bark  thick- 
ness were  not  different  from  zero. 

Thus,  for  sampling  length  of  egg  galleries  (which  are  pro- 
portional to  initial  beetle  population),  we concludethattree 
diameter  and  sample  height,  not  bark  thickness,  need  be 
considered.  As  we  discuss  later,  however,  mortality  factors 
such  as  C.  vancouverensis  are  affected  by  bark  thickness. 

In  sampling  Douglas-fir  beetles  it  is  less  costly  to  take 
samples  from  the  more  accessible  lower  part  of  standing 
trees.  Would  this  result  in  a  biased  estimate  of  the  density 
of  bark  beetles?  In  table  2  the  density  at  1.5  m  is  about  one- 
half  the  average  tree  density.  The  density  at  4.6  m  is  very 
close  to  the  average  density  of  egg  galleries,  and  thesimple 
correlation  of  the  4.6  m  sample  with  the  tree  average  is 
Ryj^  =  0.640,  (Ryx  =  0.410).  Therefore,  average  Douglas- 
fir  beetle  attack  density  per  tree  is  estimated  better 
at  the  4.6  m  height  (a  slightly  lower  height  may  suffice 
[Furniss  1962b]). 

In  order  to  relate  4.6  m  samples  to  the  tree  as  a  whole, 
however,  one  needs  a  regression  estimator  (Kish  1965). 
From  a  regression  of  the  average  beetle  density  per  tree  on 
the  density  at  the  4.6  m  level,  the  following  relationship  was 
obtained: 

Y,  =  0.750  +  0.580  X,, 

where  Vj  is  average  beetle  density  in  the  tree,  X,  is  the  beetle 
density  in  the  4.6  m  level  sample. 


Table  1  .—Length  of  egg  galleries  (dm)  and  maximum  bark  thickness  (cm)  by  sample  height 
(m)  for  5  sample  trees  in  each  of  5  diameter  classes.  The  diameter  (cm)  for  each 
tree  at  1.5  m  height  is  also  presented 

Sample 

height   Average  diameter  (cm)  at  1.5  m  height  

(m)  27.9  43.2  63.5  71.1  81.3 


32.0 

1.27(2.06) 

29.0 

1.02(0.48) 

25.9 

1.27(0.25) 

22.9 

1.02(2.46) 

1.52(1.37) 

1.78(3.53) 

19.8 

1.02(2.29) 

1.27(2.34) 

1.52(1.07) 

16.8 

1.02(2.79) 

1.78(2.64) 

1.78(0.81) 

13.7 

0.76(0.0) 

1.52(2.34) 

1.78(2.54) 

2.29(2.01) 

10.7 

1.02(1.80) 

1.27(2.41) 

2.29(1.80) 

2.03(3.00) 

7.6 

1.02(0.0)^ 

1.27(3.23) 

2.03(1.45) 

2.03(1.45) 

2.29(1.32) 

4.6 

1.02(.91) 

1.52(1.73) 

1.78(1.68) 

3.30(1.32) 

2.54(1.32) 

1.5 

1.27(1.40) 

2.79(0.0) 

3.30(0.76) 

5.08(0.0) 

5.59(1.78) 

'Maximum  bark  thickness  (cm)  and  (length  of  egg  galleries  (dm)). 
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Table  2."Average  and  variance  of  length  of  Douglas-fir  beetle 
egg  galleries/sample  per  tree  at  the  two  lowest  heights 

 I  

Height  Length  of  galleries  (dm)  Variance 

1.5  m  0.81  0.68 

4.6  m  1.60  .57 
Average  for  tree  1-67  .46 


Predators  and  Parasites 

Sampling  the  Douglas-fir  trees  at  the  4.6  m  height 
may  provide  useful  estimates  of  the  density  of  bark  beetles 
per  tree,  but  what  of  predator  and  parasite  densities? 

Table  3  shows  that  for  some  predators  and  parasites 
(including  the  two  most  abundant)  the  4.6  m  sample  is 
not  a  good  estimate  of  the  average  density  per  tree.  For  C. 
vancouverensis  the  4.6  m  sample  greatly  underestimates 
the  population  (all-tree  sample)  density.  For  Medetera  spp. 
the  4.6  m  sample  overestimates  the  all-tree  sample  by  a 
large  amount.  Roptrocerus  sp.  and  Cecidostiba  sp.  also 
have  much  larger  densities  per  tree  based  on  all  samples 
than  the  density  based  on  4.6  m  samples.  However,  these 
last  two  occur  in  very  low  frequency  and  are  probably 
mainly  parasites  of  other  (secondary)  bark  beetles  high 
in  the  tree. 

Because  C.  vancouverensis  and  Medetera  spp.  were  the 
most  abundant  parasite  and  predator,  their  distributions 
on  the  tree  will  be  evaluated  further. 


Table  3.-Average  number  of  predators  and  parasites  per  sample 
per  tree  and  per  sample  at  4.6  m  height 


Species 

Number 
per  sample 
per  tree 

Number  per 
sample  at 
4.6  m  height 

Enoclerus  sphegeus 

0.4973 

0.5833 

Thanasimus  undatulus 

.0699 

.0416 

Temnochila  chlorodia 

.3333 

.4375 

Coeloides  vancouverensis 

3.4919 

1 .7500 

Roptrocerus  sp.  and  Cecidostiba  sp. 

.4274 

.2083 

Medetera  spp. 

2.1693 

3.5000 

Total 

7.4462 

6.9375 

COELOIDES  VANCOUVERENSIS 

Of  the  372  samples  taken,  80  contained  C. 
vancouverensis.  Of  the  48  trees  sampled,  21  contained 
C.  vancouverensis.  One  hundred  and  twenty-three 
C.  vancouverensis  were  found  on  samples  in  one  tree 
and  50  occurred  on  one  sample.  This  indicates  a 
contagious  distribution.  We  herein  explore  some  reasons 
for  this  distribution. 

Coeloides  vancouverensis  deposits  eggs  singly  through 
the  bark  onto  a  2nd,  3rd,  or  4th  instar  D.  pseudotsugae 
larva  (Ryan  and  Rudinsky  1962).  Larvae  of  this  parasite 
were  found  beginning  July  12,  12  days  after  start  of 
sampling.  The  female  is  excluded  from  ovipositing  in 
those  areas  of  the  tree  where  overlying  bark  is  thicker  than 
the  5  mm  length  of  ther  ovipositor.  Density  and  distribution 
of  host  larvae,  also  govern  the  location  and  abundance  of 
C.  vancouverensis. 


To  model  the  distribution  of  C.  vancouverensis  at  a 
sample  location,  a  logistic  function  was  chosen.  The 
response  variable  (Y)  was  recorded  as  1  if  the  parasite 
was  present;  0  if  the  absent  on  a  sample.  The  logistic 
function  which  defines  the  probability  of  the  presence  of  a 
given  insect  is  given  by  V  =  P(n  >  0)  =  (1  +  EXP[-BX])-i, 
where  n  is  the  number  of  C.  vancouverensis  (an  integer 
greater  than  or  equal  to  zero).  EXP  is  the  exponential 
function.  X  is  a  vector  of  independent  variables  (bark 
thickness,  height,  etc.).  B  is  a  vector  of  coefficients  to  be 
estimated. 

Walker  and  Duncan  (1967)  developed  an  efficient 
method  of  estimation  for  the  coefficient  vector  B.  The 
algorithm  has  been  packaged  as  a  computer  program 
in  Fortran  IV  (Hamilton  1974).  Using  that  procedure,  the 
following  model  was  derived  for  individual  samples: 

=  (1  +  EXP[f(X.,8)]-i  +  ej 

where 

f(Xj,8)  =  8o  +  SiXij  X  82X2  +  83X3.  +  84X4J 

and  ej  is  the  error  associated  with  the  yth  sample. 
The  variables  used  have  the  following  identification: 

Vj     =  the  presence  of  C.  vancouverensis  (1  =  present, 

0  =  absent), 
X-,.  =  date  of  the  /th  sample, 

X2j  =  length  of  egg  galleries  (dm)  in  the  /th  sample, 
X3j  =  maximum  bark  thickness  (cm)  of  the/th  sample, 
X4j  =  number  of  ventilation  holes  in  the  /'th  sample. 

The  estimates  of  these  coefficients  based  on  this  data 
set  are: 

i)0  = -13.6380  ^3  = -0.6450 

b^  -    0.0650  b4  =  0.1680 

t>2  =  0.1299 

The  degree  of  fit  is  seen  in  table  4.  Here  the  observations 
are  grouped  into  classes  according  to  their  predicted 
values.  For  example.  Class  1  contains  only  one  observa- 
tion. Of  the  372  samples,  the  model  predicted  a  probability 
of  0.01  or  less  C.  vancouverensis  for  only  one  sample. 
This  sample  did  not  contain  C.  vancouverensis.  By 
contrast,  Class  20  contains  those  samples  for  which  the 
predicted  probability  of  containing  C.  vancouverensis  is 
0.90  <  0.95.  This  class  contains  10  samples  of  which  9 
contain  C.  vancouverensis.  The  excellent  agreement  of 
expected  and  observed  C.  vancouverensis  is  indicated  by  a 
chi-square  statistic  of  0.01  for  this  class.  An  overall  chi- 
square  for  the  21  classes  is  13.60,  which  indicates  good 
degree  of  fit  for  the  logistic  function. 
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Table  4.— Number  of  samples  with  C.  vancouverensis  as  a  function  of  the  expected  frequency 
and  chi-square  values  for  degree  of  agreement  between  the  expected  and  observed 
frequencies 


Probability  of 

Number  of  samples 

Expected  number 

presence  of 

Number  of 

with 

of  samples  with 

Chi- 

Class 

C.  vancouverensis 

samples 

C.  vancouverensis 

C.  vancouverensis 

square 

1 

0.0  - 

0.01 

1 

0 

0.0 

0.01 

2 

.01  - 

.05 

9 

0 

.3 

.27 

3 

.05  - 

.10 

23 

0 

1.7 

1.73 

4 

.10  - 

.15 

35 

"I 

4.4 

2.63 

5 

.15  - 

.20 

33 

4 

5.8 

.55 

6 

.20  - 

.25 

27 

6 

6.1 

.00 

7 

.25  - 

.30 

21 

8 

5.8 

.87 

8 

.30  - 

.35 

19 

11 

6.2 

3.83 

9 

.35  - 

.40 

14 

7 

5.2 

.59 

10 

.40  - 

.45 

9 

2 

3.8 

.88 

11 

.45  - 

.50 

8 

5 

3.8 

.38 

12 

.50  - 

.55 

5 

2 

2.6 

.15 

13 

.55  - 

.60 

18 

13 

10.3 

.70 

14 

.60  - 

.65 

15 

11 

9.4 

.29 

15 

.65  - 

.70 

20 

14 

13.5 

.02 

16 

.70  - 

.75 

18 

12 

13.0 

.09 

17 

.75  - 

.80 

30 

21 

23.2 

.23 

18 

.80  - 

.85 

34 

29 

28.0 

.03 

19 

.85  - 

.90 

18 

14 

15.7 

.20 

20 

.90  - 

.95 

10 

9 

9.2 

.01 

21 

.95  - 

1.00 

5 

4 

4.9 

.16 

Total 

13.60 

Figure  17  shows  the  relationship  of  date  and  the 
probability  of  C.  vancouverensis  being  present,  if  other 
variables  in  the  model  are  held  at  their  average  values.  This 
function  reflects  the  rapid  increase  in  frequency  of  C. 
vancouverensis  during  midsummer  (fig.  18).  The 
probability  of  presence  of  C.  vancouverensis  is  also 
positively  correlated  with  the  average  length  of  egg 
galleries  (fig.  19).  Bark  thickness  strongly  affects  the 
probability  of  occurrence  of  C.  vancouverensis  (fig.  20). 
Samples  at  4.6  m  height  contain  thicker  bark  than  samples 
taken  higher  up  the  tree,  and  this  explains  why  the  4.6  m 
sample  under  estimates  the  average  density  of  C. 
vancouverensis. 

MEDETERA  SPP. 

Medetera  spp.  were  the  most  numerous  predators  in  the 
trees  and  occurred  over  the  entire  samplmg  time  period 
(fig.  21).  They  were  found  more  frequently  atthe4. 6m  level 
than  on  the  average  of  all  samples  from  the  tree.  Greater 
moisture  in  thick  bark  portions  of  the  tree  probably  protect 
the  larvae  from  desiccation. 

Again  the  logistic  model  can  be  employed  to  determine 
where  the  Medetera  spp.  are  most  likely  to  occur.  The 
model  is  of  the  form: 


=  P(nj>0)  =  (1  +  EXP[Xj,8])-i  +  ej 

where  the  function  f(X^,B)  will  contain  both  linear  and 
quadratic  terms  and     is  the  random  error  component. 

The  following  variables  were  used  in  the  model: 

Vj  =  the  presence  of  Medetera  spp.  (1  =  present,  0  = 
absent), 


X 

X2 
X3 

X4 

^5 
^6 
Xj 

Xs 

X9 


1^  -  the  date  of  the  /'th  sample, 

2 


^1 


J' 


the  height  (m)  of  the  /'th  sample, 

x|. 

J,' 

the  maximum  bark  thickness  (cm)  of  the  /th  sample, 

the  number  of  ventilation  holes  in  the  /th  sample, 

2 

the  diameter  (cm)  of  the  tree  at  the  height  where 
the  /th  sample  was  taken. 


The  sample  estimates  of  the  coefficients  for  the  above 
model  are: 
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Figure  1 /.--Probability  of  C.  vancouverensis  on  a  sample  as  a  function  of  day  of 

the  year. 
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Figure  1 8. —Average  number  of  C.  vancouverensis  per  sample  per  tree  by  day 
of  the  year  on  which  the  tree  was  sampled. 
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Figure  20.— Probability  of  C.  vancouverensis  on  a  sample  as  a  function  of  maximunn 

bark  thickness. 
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Figure  21  .--Average  number  of  Medetera  spp.  per  sample  per  tree  by  day  of  the 
year  on  which  the  tree  was  sampled. 


bo  =  -65.0992 
b^  =  +0.57756 
b2  =  -0.001366 
b3  =  +0.169242 
b4  =  -0.0081439 


bs  =  +1.16607 
be  =  -0.183806 
bj  =  +0.620687 
bQ  =  -0.057379 
bg  -  +0.05076 


The  model  closely  fits  the  observed  distribution  of 
Medetera  (chi-square  12.82,  table  5).  In  examinating  the 
form  of  the  model,  we  will  discuss  the  sample  variables 
individually,  then  jointly. 

The  probability  of  Medetera  occurring  on  a  sample  in 
relation  to  changes  in  maximum  bark  thickness  is  shown 
in  figure  22.  All  other  variables  in  the  model  are  held  at 
their  average  values.  The  form  of  the  curve  fits  our  obser- 
vations that  Medetera  were  rare  at  upper  heights  where 
bark  is  thinner.  We  also  know  that  Medetera  were  more 
abundant  than  average  at  the  4.6  m  level.  The  peak  of  the 
curve  occurs  at  about  3.25  cm,  half  way  between  the 
average  maximum  bark  thickness  for  the  4.6  m  and  1.5  m 
samples.  The  decline  in  abundance  of  Medetera  where 
bark  was  thicker  than  3.25  cm  is  due  to  the  influence  of 
1.5  m  samples  (av.  3.90  cm  max.  bark  thickness),  there 
being  fewer  Douglas-fir  beetle  attacks  and  progeny  at 
that  sample  height  (Furniss  1962b).  Medetera  decreased 
in  density  with  sample  height  (fig.  23)  and  increased  with 
diameter  (fig.  24)  at  the  sample  location. 
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Table  5.~Number  of  samples  containing  Medetera  spp.  as  a  function  of  the  expected 
frequency  and  chi-square  values  for  degree  of  agreement  between  the 
expected  and  observed  frequencies 


Probability  of 

rMUmuer 

numocr  ui 

cxpecieu  numuer 

presence  of 

ot 

samples  with 

of  samples  with 

oni- 

Class 

Medetera  spp. 

samples 

ModotorB  spp. 

Medetera  spp. 

square 

1 

0.0  - 

0.01 

A 

4 

U 

n  n 
U.U 

n  no 

2 

.01  - 

.05 

■1  c 

1  D 

-4 

1 

.4 

1  .vd 

o 
O 

.05  - 

.10 

Q 
0 

Z 

.b 

o.yb 

A 

4 

.10  - 

.15 

1  3 
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1  .D 

O    A  CZ 

3. 4b 
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.20 

1  b 

o  c 
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.  1  o 
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.20  - 
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1  4 

d. 
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15 
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.70 

27 

18 

18.2 
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16 

.70  - 

.75 

33 

23 

23.9 

.04 

17 

.75  - 

.80 

37 

29 

28.6 

.00 

18 

.80  - 

.85 

25 

23 

20.7 

.26 

19 

.85  - 

.90 

24 

22 

20.9 

.06 

20 

.90  - 

.95 

14 

13 

12.9 

.00 

21 

.95  - 

1.00 

1 

1 

1.0 

.00 

Total 

12.82 
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Figure  22. -Probability  of  Medetera  spp.  on  a  sample  as  a  function  of  maximum  bark 

thickness. 
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Figure  23. --Probability  of  Medetera  spp.  on  a  sample  as  a  function  of 

sample  height. 
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Figure  24. —Probability  of  Medetera  spp.  on  a  sample  as  a  function  of  tree 
diameter  at  the  sample  location. 
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The  relationship  between  Medetera  and  Douglas-fir 
beetle  ventilation  holes  is  shown  in  figure  25.  Low  numbers 
of  ventilation  holes  probably  indicate  a  low  Douglas-fir 
beetle  density  and  few  opportunities  for  the  predators  to 
oviposit  and  to  locate  prey  in  order  to  survive.  On  the  other 
hand,  very  numerous  ventilation  holes  may  have  exceeded 
the  capability  of  this  particular  Medetera  population  to 
exploit  them  as  oviposition  sites. 

The  effect  of  date  of  sample  is  seen  in  figure  26.  This 
function  may  be  influenced  by  a  few  trees  sampled  in  mid- 
summer that  had  a  very  high  density  of  Medetera  spp. 
The  subsequent  decline  in  Medetera,  however,  might  be 
due  to  drying  of  phloem  as  the  season  progressed,  which 
would  likely  take  a  toll  of  Medetera  larvae,  especially 
in  areas  of  thinner  bark. 

The  joint  effect  of  the  physical  variables  that  describe 
the  sample  location  on  a  tree  are  seen  in  figure  27.  Here 
date  and  number  of  ventilation  holes  are  held  at  their 
average  values.  The  diameter,  maximum  bark  thickness, 
and  height  are  values  for  this  average  tree.  At  heights  up 
to  13.7  m,  the  probability  of  finding  Medetera  spp.  is  high 
(73-86  percent).  As  the  samples  are  taken  at  heights  above 
14  m,  the  probability  that  Medetera  spp.  are  present 
declines. 
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Figure  25.— Probability  of  Medetera  spp.  on  a  sample  by  number  of 
Douglas-fir  beetle  ventilation  holes  per  sample. 
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Figure  26.-Probabllity  of  Medetera  spp.  on  a  sample  by  day  of  the  year  on 

which  the  sample  was  taken. 
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Figure  27.— Probability  of  Medetera  spp.  on  samples  on  a  tree  using 
all  sample  variables  in  the  prediction  function. 
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Effect  of  Predators  and  Parasites  on  Douglas-fir  Beetle 
Density 

How  many  DFB  progeny  do  predators  and  parasites  kill? 

Because  progeny  are  consumed,  mortality  cannot  be 
measured  directly.  But  mortality  can  be  estimated  if  the 
number  of  Douglas-fir  beetles  per  sample  is  viewed  as  a 
function  of  length  of  egg  gallery  and  number  of  predators 
and  parasites.  This  relationship  can  be  expressed  as 
follows: 

V^j-  =  number  of  the  Douglas-fir  beetle  found  in  the  /th 
sample  of  the  /th  tree, 

X|^,j  =  number  of  the  kth  predators  or  parasites  found  in 
the  /th  sample  of  the  /th  tree,  (k  =  1,  2,.  .  .,  7), 

X31  j  =  length  of  egg  galleries  (dm)  found  in  the  /th  sample 
of  the  /th  tree. 

The  predators  and  parasites  are  identified  as  follows: 


X-||  j  =  E.  sphegeus 

X2ij  =  T.  undatulus 
X3jj  =  T.  chlorodia 
X4i,  =  C.  vancouverensis 


X5|j  =  Roptrocerus  sp.  and 
Cecidostiba  spp. 
Xeij  =  Medetera  spp. 
Xy,;  =  Belosta  albipilosa 


The  coefficients  (S's)  are  estimated  by  linear  regression. 
Regression  coefficient/significant  level  (using  the  t-test 
are): 

So  =  +4.472/.001  65  =  -0.797/.081 

St  =  +0.062/.891  Be  =  -0.274/.074 

82  =  -3.690/.047  Bi  =  -1.834/.001 

83  -  -0.670/.299  88  =  +3.901 6/.001 

84  =  -0.392/.001 

The  degree  of  fit  of  this  function  is  indicated  by  multiple 
correlation  coefficient  R2  =  0.24  and  a  standard  deviation  of 
9.325.  The  regression  function  has  accounted  for  24 
percent  of  the  variation  in  the  number  of  Douglas-fir 
beetles  per  sample.  Most  coefficients  in  the  equation  are  of 
correct  sign  as  indicated  by  their  supposed  effect,  but 
the  magnitude  of  the  coefficients  is  smaller  than  expected. 
For  example,  C.  vancouverensis  is  a  parasite,  so  the 
coefficient  84  should  approach  -1.0  (one  less  Douglas-fir 
beetle  per  sample  for  each  C.  vancouverensis  per  sample). 
The  coefficient  81  is  of  the  wrong  sign  and  is  essentially 
zero.  Both  Bi  and  S3  are  not  significantly  different  from 
zero  using  Student's  t-test  and  a  =  0.10  level  for  testing. 

If  coefficient  83  for  length  of  egg  galleries  was  closer 
to  11.26,  the  coefficients  of  the  other  terms  could  be 
interpreted  as  the  efficiency  of  that  predator  or  parasite. 
The  ratio  of  DFB  progeny  to  length  of  egg  galleries,  under 
low  pressure  from  entomophagous  insects,  is  1 1 .26.  Under 
this  condition  mortality  would  be  BqXq,^  -  Vjj.  At  any  value 
of  88,  however,  the  coefficients  still  measure  the  relative 
effectiveness  of  one  predator  to  another. 


SAMPLING  RECOMMENDATIONS 

Measurements  of  Douglas-fir  beetle  attacks  and 
survival  provide  useful  population  indices,  which  foretell 
tree  mortality  prior  to  fading  of  trees  (Furniss  and  others 
1979).  When  tempered  with  weather  data  and  stand 
susceptibility,  population  measurements  could  alert 
forest  managers  of  the  need  for  intensive  surveys  and 
management  actions  up  to  a  year  ahead  of  present  aerial 
detection  surveys. 

But  where  should  samples  be  located?  When  should 
sampling  occur  and  how  many  samples  are  needed?  Some 
of  these  questions  have  been  addressed  before  (Furniss 
1962a,  b;  1964).  Samples  taken  closest  to  the  time  of 
emergence  of  beetle  progeny  will  be  correlated  best  with 
surviving  population  level  (Furniss  1964).  Fall  sampling, 
however,  is  a  practical  compromise  because  mortality 
rate  is  fairly  stable  by  then  and  work  is  more  easily 
accomplished  because  roads  are  more  accessible  and 
weather  is  better  than  in  early  spring. 

There  is  no  single  location  in  a  tree  that  is  representative 
of  beetle  attacks,  progeny,  and  entomophagous  insects. 
Due  to  the  particular  distributions  of  the  Douglas-fir 
beetle  and  its  two  most  abundant  parasites  and  predators, 
future  sampling  should  consist  of  a  minimum  of  three 
locations  per  tree:  (1)  4.6  m,  (2)  midpoint  of  the  infested 
stem,  and  (3)  about  3  m  below  the  top  of  the  infested  stem 
to  avoid  the  characteristic  resinous,  unsuccessful  zone 
occurring  near  the  top  (Furniss  1962b).  Using  those 
sample  locations,  we  calculate  that  it  would  require  4 
and  30  trees  to  estimate  average  densities  of  progeny 
and  egg  galleries,  respectively,  per  sample  per  tree,  with 
CV-    =  0.20. 

If  trees  cannot  be  felled  for  sampling,  an  alternative  might 
be  to  sample  at  4.6  m  height  with  a  ladder  and  adjust 
average  beetle  density  with  a  regression  estimator,  which 
in  our  study  was  /|  =  0.750  +  0.580  X,.  In  that  case,  we 
would  have  needed  to  sample  6  and  30  trees  to  estimate 
progeny  and  egg  galleries,  respectively,  at  the  4.6  m  height 
with  a  CVj^  =  0.20.  The  sampler  should  note,  however,  that 
the  regression  estimator  may  differ  in  other  infestations 
requiring  additional  sampling  to  calculate  the  regression 
coefficients. 

The  question  of  sampling  individual  species  of 
entomophagous  insects  is  more  difficult.  The  problem  is 
complicated  by  low  densities  of  some;  seasonal  variations 
in  their  presence  or  absence;  and  their  distributions. 
Accordingly,  the  sampler  cannot  assess  them  collectively 
at  any  one  time  or  location.  In  order  to  measure  their 
effects  (density  of  surviving  beetle  progeny),  however, 
sampling  should  occur  after  mortality  has  leveled  off  at 
the  end  of  summrer. 
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Nonetheless,  if  we  take  our  samples  from  the  three 
recommended  heights,  the  numbers  of  trees  required  for 
a  CV^^  =  0.20  for  the  seven  categories  of  entomophagous 
insects  are: 


Coeloides  47  (37)^ 

Medetera  28  (27) 

Ecoclerus  67  (64) 

Thanasimus  595  (132) 

Temnochila  37  (29) 

Roptrocerus,  Cecidostiba  85  (76) 

Belosta  67  (56) 


Tor  comparison,  numbers  in  parentheses  are  samples  needed  using 
data  from  3-m  Intervals  instead  of  three  samples  per  tree. 

CONCLUSIONS 

The  mathematical  model  based  on  the  logistic  function 
represents  realistically  the  seasonal  abundance  of 
Douglas-fir  beetle  progeny  in  relation  to  regulatory  factors 
involving  entomophagous  insects  and  sample  charac- 
teristics. The  model  is  useful  for  demonstrating  the 
general  importance  of  parasites  and  predators  in  the  post- 
egg-hatch  period.  For  example,  average  mortality  of 
Douglas-fir  beetle  progeny  after  egg  hatch  in  successfully 
attacked  trees  was  58  percent.  Further  improvement  in 
the  model  should  be  directed  toward  separation  of  the 
proportional  mortality  due  to  entomophagous  insects  and 
that  due  to  competition  among  progeny. 

The  information  presented  on  the  temporal  and 
spatial  distribution  in  trees  of  the  Douglas-fir  beetle 
and  its  major  enemies  will  facilitate  evaluation  of  in- 
festations and  development  of  an  improved  pest  manage- 
ment program.  Such  a  program  will  need  to  include  infor- 
mation on  susceptibility  of  surrounding  trees  and  stands 
because  that  factor  limits  the  amount  of  mortality  that  a 
particular  population  can  cause.  Research  at  this  Station, 
in  cooperation  with  Forest  Service  Region  1,  and  Idaho 
Department  of  Lands,  is  being  directed  to  meet  that  need. 
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The  Intermountain  Station,  headquartered  in  Ogden, 
Utah,  is  one  of  eight  regional  experiment  stations  charged 
with  providing  scientific  knowledge  to  help  resource 
managers  meet  human  needs  and  protect  forest  and  range 
ecosystems. 

The  Intermountain  Station  includes  the  States  of 
Montana,  Idaho,  Utah,  Nevada,  and  western  Wyoming. 
About  273  million  acres,  or  85  percent,  of  the  land  area  in  the 
Station  territory  are  classified  as  forest  and  rangeland.  These 
lands  include  grasslands,  deserts,  shrublands,  alpine  areas, 
and  well-stocked'  forests.  They  supply  fiber  for  forest  in- 
dustries; minerals  for  energy  and  industrial  development;  and 
water  for  domestic  and  industrial  consumption.  They  also 
provide  recreation  opportunities  for  millions  of  visitors  each 
year. 

Field  programs  and  research  work  units  of  the  Station 
are  maintained  in: 

Boise,  Idaho 

Bozeman,  Montana  (in  cooperation  with  Montana 
State  University) 

Logan,  Utah  (in  cooperation  with  Utah  State 
University) 

Missoula,   Montana   (in   cooperation   with  the 
University  of  Montana) 

Moscow,  Idaho  (in  cooperation  with  the  Univer- 
sity of  Idaho) 

Provo,  Utah  (in  cooperation  with  Brigham  Young 
University) 


Reno,  Nevada  (in  cooperation  with  the  University 
of  Nevada) 


